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In the 2-diinensional anisotropic classical Heisenberg model with Xy-symmetry there are nonpla- 
nar vortices which exhibit a localized structure of the 2-components of the spins around the vortex 
center. We study how thermal noise induces a transition of this structure from one polarization to 
the opposite one. We describe the vortex core by a discrete Hamiltonian and consider a station- 
ary solution of the Fokker-Planck equation. We find a bimodal distribution function and calculate 

,, the transition rate using Langer's instanton theory (1969). The result is compared with Langevin 

**! I dynamics simulations for the full many-spin model. 

Tij- ■ PACS: 05.40.-hj, 75.10.Hk, 02.50.Ey, 75.30.Gw 
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I. INTRODUCTION 

There are several classes of quasi-2D (two-dimensional) magnetic materials for which the ratio of inter- to intra- 
plane magnetic coupling constants is typically 10^'^ — 10^^: (1) Layered magnets [|l|,||,|[Q], like K2CUF4, Rb2CrCl4, 
(CH3NH3)2CuCl4, and BaM2(X04) with M=Co, Ni, ... and X=As, P, ... (2) C0CI2 graphite intercalation compounds 
[g[, (3) magnetic lipid layers, like manganese stearate Q. Many of these materials can be described by the classical 
2D Heisenberg model with XY- or "easy-plane" symmetry (section II). 

In this model vortices play the decisive role: they are responsible for a topological phase transition [^|j|] at the 
Kosterlitz-Thouless temperature Tc and, above Tc, for "central peaks" in the dynamic form factors for the spin 
correlations. The central peaks were observed in inelastic neutron scattering experiments pJlQ,0,n3,n3| and in 
combined Monte Carlo/Spin Dynamics simulations [IJ jlj,|l^,|l^,|l8|,^ . The observed central peaks agree qualitatively, 
partially even quantitatively, with the central peaks which were obtained by a vortex-gas approach [ |l^ , p^ , p^|Jl^Jl^ , [l9| . 

There are two types of static vortex solutions whose structure and energy differ, depending on the anisotropy of the 
Heisenberg exchange interaction [fiw. For strong anisotropy (i. e., if the anisotropy parameter S exceeds a threshold 
6c) only planar vortices are stable for which all spins are lying in the easy plane (xy-plane). For weak anisotropy 
(0 < (5 < Sc) only nonplanar vortices are stable, which exhibit a localized structure of the z-components of the spins 
around the vortex center. In addition to the vorticity q = ±1, ±2, ..., the nonplanar vortices have a second topological 
charge p. It is denoted "polarization" because its sign determines the side of the xy— plane to which the out-of-plane 
vortex structure points. The planar vortices can be considered as having p = 0. 

The product qp of the topological charges determines the dynamics because the vortices are subject to a "gyrocou- 
pling force" G x V, which is formally equivalent to the Lorentz force |2y,|2^: V is the velocity of the vortex center, 
but instead of an external magnetic field we have here an intrinsic quantity, produced by the vortex itself and carried 
along with it: The "gyrovector" G = 2TTqpCz which is orthogonal to the xy-plane. The formula for G was derived in 
the continuum limit and, strictly speaking, G is conserved only in this limit. Nevertheless, spin dynamics simulations 
for 1 or 2 vortices showed that the direction of G (or the sign of p, because q is always conserved) does not change 
during the simulation p^ , p3t . 

However, we know so far three exceptions, i. e. situations in which the out-of-plane vortex structure can suddenly 
make a transition from one polarization to the opposite one. As the direction of G is reversed, this has a drastic 
effect on the dynamics: The direction of the gyrocoupling force is also reversed, which means that the direction of 
the vortex motion is reversed, too. The three transition mechanisms are: 

(1) Interaction with spin waves. The easiest way to see this is to use "dirty" initial conditions for the spin dynamics 
simulation [p4|: E. g., a structure which is not a good approximation to the 1-vortex solution (this solution can be 
obtained numerically by an iteration procedure |24| ). Then many spin waves are radiated at the beginning of the 
simulation, while the approximate vortex structure adapts to the lattice and becomes a "good" solution (numerically 
identical to the above solution obtained by iterations). The emitted spin waves form a magnon gas; i. e. the vortex 
moves in a kind of magnon thermostat and transitions to the opposite polarization occur with a certain probability 
which depends on how dirty the initial condition was. 

(2) An ac magnetic field. If the amplitude of a field which rotates in the easy plane is larger than a threshold value, 
a transition to the opposite polarization occurs. In contrast to (1), the reverse process does not occur because the 
field breaks the symmetry of the two polarizations. This will be the subject of a forthcoming paper. 

(3) Thermal noise. This has some similarity with (1), although that is a deterministic zero-temperature effect. In 
section II we implement white noise into the microscopic equations (the Landau-Lifshitz eq.) by adding stochastic 
magnetic fields to the local fields in which every spin precesses. In this way we model the interactions of the spin 
degrees of freedom with thermostat degrees of freedom (magnons, phonons etc.). We consider a stationary solution Vst 
of the Fokker-Planck equation, using a reduced Hamiltonian which models the vortex core. Such a core Hamiltonian 
was used in [^,^ for the calculation of 6c- For a certain parameter range, Vst exhibits two maxima (for the two 
possible polarizations of a nonplanar vortex) and a saddle point (corresponding to the planar vortex). 

In section HI we calculate the probability flux over the region around the saddle point using Langer's instanton 
theory p^ . Here we use the fact that for S ^ Sc there is a soft mode among the normal modes which were obtained 
numerically for a system with one vortex [Eq|. 

Finally, our prediction for the transition rate is tested by Langevin dynamics simulations, i. e. by integration 
of the stochastic Landau-Lifshitz equation. For these tests the design of the simulations, including the choice of the 
parameter ranges, turns out to be decisive. 



II. HAMILTONIAN AND THE FOKKER-PLANCK EQUATION 

We consider a Heisenberg model with XY- or easy-plane symmetry with classical spins Sn located on the sites 
n = {ux, Uy) of a square lattice 

H--Y.JK {Sn Sl_^ + SI Sl_^ + A S^ 5^_^) (1) 

n,A 

where A is the anisotropy parameter (0 < A < 1), J^= J is the exchange integral and A — (A^, Aj,) is a vector 
which connects a spin with its nearest neighbors (Aa; = ±l,Aj, = or A^ — ztljA^; = 0). The spin dynamics is 
governed by the Landau-Lifshitz equation. Since we want to study the interaction with thermal noise, we implement 
a noise and a damping term 

Here we have added a stochastic magnetic field /i,-j(i) to the local field -^S- in which the spin Sn precesses. Since hfi 

is multiplied with Sfi, this means multiplicative noise. 

Another way to obtain the stochastic term in Eq. ^ consists in adding to the Hamiltonian interactions between 
the spins and local stochastic magnetic fields 

Vit) = -Y,hn{t)Sn- (3) 

n 

We use Gaussian white noise with 

<hm> =0 , 

< /i^(t) h^lit') > =2D^ 6ftn' Sec' Sit - t') , (4) 

where D^ is the variance of the noise. In order to preserve the isotropy in the easy plane we demand 

D,^Dy = D . (5) 

The last term in Eq. (0) represents damping in the Landau-Lifshitz form (see m). An alternative would be the 
Gilbert damping which yields the same results, however, as we will use only very small damping coefhcients. 

It is convenient to use a representation for the classical spin vector Sfi in terms of two angles of rotation 6fi and $s 

Sfi — S{sin 9,1 cos $,-5, sin 9,^ sin <i>,-j, cos 9fi} (6) 

The variables Mfi — cos 9f{ and <I>s constitute a pair of canonically conjugated variables, which means that in the 
no-damping case (7 = 0) 



d^n _ d{H + V) 

1 

(7) 



dt dMn 

dMn d{H + V) 



dt d^fi 

Here 

H = -JY^ {xMnM._^ + PnP,^_^ cos($« - $-_^; 
^(^) = - E ( ''^nit) Mft + Pn (hUt) cos($s) + /i|(f) sin($s))) (8) 



is the Hamiltonian of the system in terms of the new variables and P,i — \/l — M|. 

Using the variables Mfj and <I>s the Landau-Lifshitz equation (0) can be written as a set of coupled stochastic 
equations 



d<^f{ _ OH 7 dH 

^ = -^-7(l-M|)^+,.(Af.,$.,.) . (9) 

where 

/,, (Ms,a>s,t) = -/i^(t) + ^ (/i|(i) cos($s) + /i^W sm($s)) , 

ffs (Ms,$s,t) = -Ps (/if^(i) sin($s) - hlit) cos($s)) (10) 

are multiplicative stochastic forces. From Eqs (Q), (|5[) and (|lO| ) we obtain 

< U {Mr, $s, t) U {Mff , ^'n' X)> =2 6{t - t') Snn' (d, + D ^^^^ cos($« - $V 

< gn {Mft, $ff , i) fffj' (A4 , ^J-j, , t') > - 2 i? 5(i - t') Snn' PnPk cos($.s - ^'n) 

< fn{Mn,^n,t) gn' {Mt,,'^L,,t') > = 2 D 6(1 - t') Snn' Mr^ smi<Pn - <i>'n) (11) 

We have introduced the stochastic magnetic fields hfi{t) to model the interaction of the spin degrees of freedom 
with thermostat degrees of freedom: phonons, electrons, other magnetic excitations etc. However, it is clear that the 
thermostat excitations are characterized by finite correlation times and a more appropriate modelling of the influence 
of these excitations would be to use colored noise. The problem under consideration is very complicated, however, 
if considered in the framework of a nonwhite-noise approach. Therefore we restrict ourselves to the case ^, where 
we understand the white noise approach as a limiting case of the colored-noise process and therefore we consider Eqs 
®' (0) ^^ Stratonovich stochastic differential equations ^. 

From Eqs (|y)-(|l^) we obtain that the equation for the probability density function 

Vimn, 4>n, t) = {Y{5 {mn - Mn{t)) S {4>n - $3(0)) (12) 

n 

has the form 

n 

As was mentioned above, the stochastic magnetic fields hfi{t) model the interaction with thermostat degrees of freedom. 
Therefore it is quite natural to demand that Eq. (|l^) has a stationary solution in the form of the Gibbs distribution 

Vst ~ exp (-^\ . (14) 



It is seen from Eq. ([T^ ) that the function (|l4| ) is a steady-state solution of the Fokkcr-Planck equation (13) when the 
fluctuation-dissipation condition 

D.^D^jT (15) 

is fulfilled. Here T is the temperature of the crystal. Under the condition (|l5| ) the Fokker-Planck equation for the 
function V has the form 

d ^-^ d { , dH 7 dH^ 1 d 

dt ^ d(j)ft \ dm,i l-mid(j)ft l-mtd(j)n 



Let us consider first the equilibrium properties of the system. We assume that a vortex is situated in the center of 
a unit cell at the origin of a coordinate system. The static in-plane vortex (tti^ = 0) is characterized by the angles 
$^ which satisfy the equation 

^sin(ci>^-<l>«_^) = 0. (17) 

A 

The $^ are approximately given by the usual in-plane vortex structure 

$^ = g arctanf ^ I (18) 

\nxj 

where rij., Uy = {2n + l)/2, n = 0, ±1, ±2, .. (the lattice constant is chosen equal to 1) and the integer q is the vorticity. 
It is known [g5| , [2^ that the in-plane vortex is stable for < A < Ac where the critical value Ac of the anisotropy 
parameter depends on the type of the lattice (e.g. for square lattices Ac — 0.703 f^^l). For A > Ac the in-plane vortex 
becomes unstable and an out-of-plane vortex is created. To gain insight as to how the temperature influences the 
stability conditions we need a reduced form of the Hamiltonian (g|) which effectively takes into account both types of 
vortices: in-plane and out-of-plane. Such an effective Hamiltonian was proposed in pq] . It was shown in Eq] that the 
dynamics of the vortex instability can be understood under the following assumptions: 

i) The in-plane angles ^^ for static in- plane and out-of-plane vortices arc given by Eq. dl7). 

ii) The deviations tp.^ = $,t — $^ of the in-plane angles from their static values are radially symmetric. They 
strongly decay with the distance r,-j — •\/(riJ^~T72p~r"(n^^~~T72)^ from the vortex center: 



■0n 



V'l, for n = ±(1/2, 1/2), ±(-1/2, 1/2) 

^2, for n = ±(3/2, 1/2), ±(-1/2, 3/2) 

^3, for n = ±(1/2, 3/2), ±(-3/2, 1/2) ^^''> 

, otherwise 



iii) The deviations of the out-of-plane components arc also radially symmetric and decay strongly 

mi, for n = ±(1/2, 1/2), ±(-1/2, 1/2) 

7712 , for n = ±(3/2, 1/2), ±(-1/2, 3/2) 

mg, for n = ±(1/2, 3/2), ±(-3/2, 1/2) ^^"^ 

, otherwise 

Under these assumptions the dynamics of the vortex core is described by the following Hamiltonian 

He ~ -4:J {A [ml + TOi (to2 + "T-a) + ™2 "t-s) 

+ COS^iPi (p2 COs(f/'l - -02) +P3 COs(f/'l - V's)) 
+ (cOS(5i + COS (52) (P2 COSV'2 +P3 COSIps) 

+ P2P3 sin(2(5i) cos(V'2 - "^s)} (21) 

with 



Pa = v/l-"^a, a = 1,2, 3 (22) 

and 

^1 = ^1/2,1/2 - *3/2,l/2 ' "^2 = ^3/2,1/2 ^ ^5/2,1/2 ' (23) 

Using the approximation ip2 — ■(/'3,"Z2 — m^ and the static in-plane angle distribution (RW, we get cos(5i = 
2/\/5, cos (52 = 8/\/65 and in this case the Hamiltonian ( |2l| ) coincides with the Hamiltonian in ||25|. 

Being interested in the distribution of the out-of-plane components ma we integrate the function (n^ with respect 
to the in-plane angles ipa- We obtain a reduced stationary probability density 

rstimi,m2,m3) = J_ g-^/'H+^i (™2+™3)+™2 ms) 

2ir 



'0 e'^ '^'"(251 )p, PS COS c^j^ (^ ^Qg g^ p^ ^p2 ^p2^ 2p2P3 COS </>) X 

/o(/3(cos(5i +cos(52)Ypi+pi + 2p2P3Cos(/)) , (24) 



where (3 — 4 J/T is a dimensionless inverse temperature. Af is the normahzation factor and Io{x) is a modified Bessel 
function ^^. 

The analysis of the function (E4) shows that it has a unique maximum at mi — rn2 — Jn^ = if the anisotropy 
parameter A is below a temperature dependent threshold value Ac(T). This case corresponds to the stable in-plane 
vortex. If 

A > \c{T) (25) 

the function (|24| ) has two maxima at mi = ± rn?, m2 — ± m2, m^ = ± 771,3 ^^'^ ^ saddle point at m,i = 7712 = ma = 0. 
In this case the probability density function p4) describes a bistable system of two out-of-plane vortex structures 
with opposite polarizations. 

Let us illustrate this statement by a crude approach when only two degrees of freedom mi and ipi are included. In 
this case the core Hamiltonian ( pi|) simplifies to 

i7c = -4 J{Am^ + 2cos(5ipiCOsV'i} (26) 

The corresponding stationary distribution Vstifniiipi) is plotted in Fig. la (without normalization). The reduced 
stationary probablity density ( p4|) has the form 

Vst{mi) = exp {I3\mj) /„ (2/3cos(5ipi) /U , 
1 

7V= / dmi exp(/3A?7i^) /q (2/3cos5ipi) (27) 

-1 

The function Vstifni) describes a bimodal distribution (see Fig. 1(b)) in the range 

a 2 r X r /i(2/3cOS(5i) 

/3 cos^(5i > A > cos (51-^-^ -{ 28 

7o(2/3cosOi) 



A more accurate approach is based on the expansion of the Hamiltonian (El^) into a series with respect to {'0}- Then in 
the harmonic approximation with respect to the ijja the stationary probability density is determined by the expression 



7'st(TOi, 7712,7773) = 

1 exp(/3iJ,({777},M-0)) 



A/" \/piP2P3(P2+P3)(sin(2(5i) (p2 + Pz) +PlC0s5i + cos 81 +C0S(S2) 

The function (|29|) describes a bimodal distribution if 



(29) 



4/3 cos ^1 (2 cos (5i + cos (^2 + 2 sin 2^1) — 2 sin 2(5i — cos (52 — 3 cos (5i 
^ ^^^ = 4/3(2 sin 2^1+ cos (52 + 2 cos (5i) ' ^ ' 

Thus the function (e3) has two maxima at 7?7i = zLm'i,m,2 = ±7772,7773 = ±7773 and a saddle point at 7771 = 7772 = 
7773 = 0. The phase diagram (the bifurcation curve A(/3)) is shown in Fig. 2. It is worth noting that for a given 
anisotropy parameter A the phase which corresponds to the in-plane vortex is always the low-tcmpcraturc phase. 

III. SWITCHING RATE 

Following Langer |g^ (see also Ref. [g9|) it is convenient to introduce a new set of variables {t^} = (771, ...ri2N) which 
consists of N out-of-plane spin deviations (771, ...?2iv) = {Tn-n} and N canonically conjugated variables {ij^^i, ■■■ri2N) = 
{(pft} and to write the Fokker-Planck equation (y^ in the forin 

where 

M,,j = TAj ~ A,j (32) 



with 

r^ = I ^v, ':3\ '.::;: ;. , , (33) 

and Aij is the following antisymmetric matrix: 



7(1 - Tyf) , for i > TV + 1 



A:j = { ~5^^3+N , J < N (34) 

, otherwise 

E{{ri}) is the Hamiltonian of the system expressed in terms of the variables {77}. 

We are interested in a switching process between the vortex states with different polarization. Therefore we 
consider the anisotropy-temperature region (see Fig. 2) where the out-of-plane vortices are stable. In this case the 
energy function £{{7]}) has a locally stable state at {rjo} (an out-of-plane vortex with positive polarization) which is 
separated by an energy barrier from another stable state {— ??o} (sm out-of-plane vortex with negative polarization). 
We assume that the system is initially prepared in a vortex state with, say, positive polarization, and we consider the 
relaxation process as an escape process from the potential well which corresponds to the vortex {tjq} neglecting the 
backward process. Another possibility to make the vortices with different polarization non-equivalent is to apply a 
constant magnetic field oriented along the hard-axis (perpendicular to the easy-plane). The in- plane vortex {fj} with 
the same vorticity as the out-of-plane vortex corresponds to the energy barrier which must be overcome. The point 
[f]} is a saddle point of E{{ri}) . 

We consider a temperature which is much smaller than the energy difference between in-plane and out-of-plane 
vortices. After having been initially in the state {rjo}, the system reaches first a quasi-equilibrium state near the 
metastable point {rjo} with the probability density V given by the Gibbs distribution 

-p ^ e^^^ with {77} ~ {ryo}. (35) 

The probability flux over the barrier is concentrated in a narrow region around the saddle point {fj} p9| ]. To obtain 
the flux, let Dni, {n, i = 1, ■.2N) be the eigenvectors of the Hessian matrix 

£^A{v})^^^E/^r^,^l^J, i,j = l,..2N (36) 

evaluated at {77} = {fj} 

2N 

Y,£^A{v})Dl,^^ilDu (37) 

and 111 are the eigenvalues. Thus the energy of the system in the immediate neighborhood of the saddle point {fj} 
can be written as 

2N 

E = E{{fj}) + -J2i,i^f (38) 



2 

1=1 



where the new variables 



2N 



6 = ^ Du iv^ - fj^) (39) 

1=1 

are the principal axes coordinates. 

Coming back to the original variables mf{ and 0fj we can say that the Hamiltonian of the system in the close vicinity 
to the in-plane vortex state can be written as 

^ -L^in — plane 

+ IjJ2{1 ^"^('^"^ - ^s+a) (^" - '/'s+a)' + cos(<I',^ - $°+^) ml - A rrm m.^^^ (40) 

ft, A 



where ipfi = <I>,-j — $^ are small deviations of the in-plane angles from their static values $^. The out-of-plane spin 
deviations mft are also assumed to be small. In this case the eigenvalues fxi correspond to the linear spin-wave 
spectrum of the system in the presence of an in-plane vortex. The normal modes were investigated in pq ] and it was 
found out that there is a particular soft mode (its frequency goes to zero for A — *■ Ac, A < Ac) which is responsible 
for the crossover from the in-p lane to the out-of-plane vortex structure. In the interval A > Ac this mode becomes 
unstable. In terms of Eq. (37) it means that the corresponding eigenvalue, say /ii, is negative. 

According to [^ (see also |2^) the rate constant for an escape from the metastable point {770} via the saddle point 
{fj} has the form 






det 



{£i{rjo})/V2^ 



T 



\ det (£({^})/V2^) 



exp 



E{{n})-E{{r^,}) \ 



(41) 



where |^| is the deterministic growth rate of the unstable mode at the saddle point. The quantity u is the negative 
eigenvalue of the following eigenvalue problem 



where Ui are eigenvectors and 



2N 

i'=i 



2N 

Mu,= ^ Di,^M,,jDi, 



(42) 



(43) 



Taking into account Eq. (|37| ) we can rewrite Eq. 

2Ar 



^) in the form 






(44) 



where Vi = '}2i=i DiaUi. 

Coming back to the original variables f/'n and m^ we obtain from Eqs. (32), (pSJ), (p4[), (|3q ) that the switching rate 
between out-of-plane vortices with opposite polarization is determined by the expression 



27r 



exp 



F 



F^. 



T 



where 



Font = Eout + T ^ In (w„(out)/T) 

m 

is the free energy of the out-plane vortex and uj.yy-^{oMi) is the m-th normal mode of the vortex. 

F^n = E,^ + T ^ In {uj^{in)/T) + T In {\u:i{in)\/T) 



(45) 



(46) 



(47) 



is an effective free energy of the in-plane vortex. In Eq. (E^) the prime means the summation over the stable modes of 
the in-plane vortex and |a;i(m)| is the modulus of the purely imaginary frequency which corresponds to the unstable 
mode of the in-plane vortex. 

The deterministic growth rate u is the negative eigenvalue of the eigenvalue problem 



E 

n' 

E 






„(i) 



\ I n' n' / 



(2)n 



drufsdrnfi' 



mii=0,*'i 



(74'^+4V) 



(1) 



.,(2) 



(48) 



where wl ' , wl, are the components of the eigenvector and we took into account that in the vicinity of the saddle point 
one can neglect the dependence on m^j in the damping constants T^. 

Let us evaluate these formulae in the crude approach already used in section 1. We consider the core dynamics taking 
into account only one pair of canonically conjugated variables mi and ipi and putting r7i2 = m^ = 0, ip2 = ips = in 
Eq. (pi|). In this case the eigenvalue problem (48) reduces to 



I'v'^^^ + 8J(X ~ Ac) {jv'^^'> + w^^^) - (49) 

where Ac = cos Si . The deterministic growth rate i/ takes on the form 

'y = 8J ({Xc - ^)7 - \V>^^1^ + 4(A - Ac)Ac') . (50) 



The out-of-plane vortex exists for A > Ac and the static value of the out-of-plane spin deviation is rrii ~ iy 1 ~ -j^- 
The Hessian matrix (pq ) is evaluated at the metastable point (out-of-plane vortex) and at the saddle point (in-plane 
vortex) which yields 

"(T|)'"(""o"'i) 



respectively. Inserting Eqs. rtSG), (|5l|) into (|4l|) yields 



£,„ -E^ 



^^,^vav+4Ac(a--a;)+7(a-2Ac) /ata:^__^ ^^^^ 



TT 

We see that in the low-damping limit the switching rate reduces to 

2J 



K^— VA2-Age- •"'^-°" (53) 

n 



while in the overdamped limit 



2Jj /A + Ac _£hl^£o^ ,_, 

K= ^(A-Ac)W— ; e T . (54) 

TT V Ac 



We note that the expressions (|5^), (^3|) and ( p4| ) are valid only when ^'" jF"^^ ':$> 1. This condition is not fulfilled 
when A ^ Ac- 

IV. LANGEVIN DYNAMICS SIMULATIONS 

In order to test our theory we have numerically integrated the stochastic Landau-Lifshitz equation (||) for a large 
square lattice in which we cut out a circle with radius L using free boundary conditions. As initial spin configuration 
we take an out-of-plane vortex with center at a distance Rq from the middle of the circle. Since the anisotropy 
parameter A should not be chosen close to Ac (see section II), the diameter 2r„ of the out-of-plane vortex structure in 
any case is considerably larger than the lattice constant. This has the advantage that the vortex can move smoothly 
over the Peierls-Navarro potential of the lattice; indeed discreteness effects are hardly visible in the motion. 

Without noise and damping the trajectory X{t) of the vortex center would be a circle with radius Rq in a first 
approximation which is given by the Thiele equation pO| 

Gx X ^ F . (55) 

The driving force F is the 2D Coulomb force between the vortex and an image vortex which is located at the distance 
L^/Rq from the circle center Je3|. The image has opposite vorticity but the same polarization as the vortex (for free 
boundary conditions). Eq. (|55|) was derived from the Landau-Lifshitz equation in the continuum limit, assuming a 
rigid vortex shape. In a better approximation the trajectories turn out to be a superposition of cycloids around the 
circular motion [p2| , but this fact seems to be unimportant for the switching process which we discuss here. 



When the damping term in (g) is included, the vortex moves outwards on a spiral Q|, until it finally reaches the 
boundary where an annihilation together with the image takes place. However, we choose an initial position far away 
from the boundary and a very small damping parameter; therefore we have plenty of time to observe the motion of 
the vortex before it gets close to the boundary. 

When the stochastic fields hnit) in (^ are included, the vortex trajectories naturally become noisy. In this case the 
variances (Xf) — (Xi)^ can be computed as a function of time and can be compared with a collective variable theory 
for finite temperature Q|, |^, |3^. This yields an effective vortex diffusion constant Dy. 

In contrast to the vorticity q, the polarization p of the vortex is not a constant of motion in a discrete system: 
The out-of-plane vortex structure can flip to the other polarization due to the stochastic fields. Then the direction of 
G = 2TTqpez is reversed and thus the direction of the vortex motion is reversed, too, as can be seen from (pa). 

In order to measure the transition rate k in the simulations it is necessary to choose carefully the parameter ranges: 
A has already been discussed above, we take A = 0.9 which is sufficiently far away from both Ac — 0.70 and the 
isotropic limit A = 1. For our circular system we choose a radius L = 24 which provides enough space for the vortex 
(the out-of-plane vortex structure should not contact the boundary even during long integration times). For the same 
reason the initial distance Rq of the vortex center from the middle of the circle should not be too large. On the other 
hand Rq should not be too small, otherwise the driving force F would not be strong enough to overcome the pinning 
forces of the lattice. Choosing Rq ~ 10 both conditions can be fulfilled, if the damping 7 is small enough. (The 
larger 7 is, the sooner the vortex reaches the boundary) . On the other hand, a small 7 means a long saturation time 
(after the start of the simulation the energy rises and saturates at a value independent of 7). For 7 > 0.002 we get 
acceptable saturation times < 300 (in units of Ti/{JS)). 

The most important parameter naturally is the temperature: For T <C -Em ~ Eout the transition rate n in (p3) is 
extremely small and thus the integration times would be much too long, which are needed to get a sufficient number 
of transition events. 

On the other hand T should not be too large, otherwise vortex-antivortex pairs appear spontaneously in the vicinity 
of the vortex. This definitely changes the translational motion of the vortex, and it is possible that the transition to 
the other polarization is influenced, too. The difference Ein — Eout can be estimated by comparing the total energies 
of our system with L = 24 in the presence of a static in-plane or out-of-plane vortex at the center of a lattice cell: 
Ein — Eout = 109.40 — 108.49 = 0.91 (in units of J). The factor in front of the exponential in Eq. ( |53| ) is approximately 
0.12 therefore 0.1 <T< 0.3 is expected to be an appropriate temperature range (The Kosterlitz-Thoulcss transition 
temperature is about 0.8 for A = 0). 

The initial spin configuration for our simulations stems from an iterative program |24| which produces a discrete 
vortex structure on the lattice (In this way we avoid the radiation of spin waves which would occur during the first 
time units if a continuum approach for the vortex structure were used). As we interprete the Landau-Lifshitz Eq. (0) 
as a Stratonovich stochastic equation and as we use multiplicative noise, we take the Heun integration scheme which 
was developed for this situation |^ , [|8[ . The spin length S is conserved in Eq. (||) and can be used as a test of the 
program, the time step is 0.01, in units of h/{JS). 

For reasons to discuss below, we have performed two different types of simulations: In type I a complete simulation 
for one temperature consists of many runs with different sequences of random numbers which produce the white noise. 
The total integration time is divided into a first part of length to (denoted as pre-run) and a second part of length t 
(denoted as main run). We choose to in the order of 1000 which is larger than the saturation time and large enough 
that the vortex has no memory of the configuration from which it started; i. e. in every run we have at the time to 
a different initial condition for the main run. Only the main runs are used for the thermal average: the average time 
T, after which the first transition of the vortex to the opposite polarization occurs, is obtained from 

N{t) = iVoe-*/^ . (56) 

Here A^o and N are the number of runs in which the vortex has made no transition until to and to + 1, respectively, r 
must be compared with the inverse transition rate k~^ — Tth from Eq. (p3|), because we work with a small damping 
parameter 7 = 0.002 (Table ||). The agreement is rather good, taking into account that we used a very crude model 
for the vortex core formed from only the four innermost spins. 

We counted only the first transitions because in our theory we have calculated the escape rate from a metastable 
state. After the first transition the vortex is typically in a different dynamical state than before, thus the probability 
for the next transition is expected to be different, too. In fact, we obtained a total number of 870 transitions in 158 
runs with t = 4000 for T = 0.15; this means that the average transition time is 917, which is about four times smaller 
than the first-transition time 4286 in Table 0. 

In a type-II simulation we only make one pre-run of length to, i. e. the main runs all start from the same initial 
condition. By taking different lengths to we can see whether t depends on to and/or the initial condition. We 



performed this type of simulations because we had some hints from the investigation of the variances |35| that a 
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certain vortex mode might be gradually excited thermally which could trigger the transition. The frequency of this 
mode is very low, namely Aw = lui — (jJ2, where u;i^2 are the eigenfrequencies of two quasi-local modes of the circular 
system with one vortex [|9J. uji,2 are identical to the frequencies of the cycloidal oscillations of the vortex trajectory 
around the mean path (see above). 

However, our type-II simulations in Table ffl do not reveal a correlation between the length t^ of the prc-run and 
the first transition time r. Nevertheless the values of r differ considerably for the different simulations. Thus we 
conclude that r depends strongly on the initial condition, which is identical for all main runs of one simulation. This 
conclusion is confirmed by looking at the first 200 time units imediately after the beginning of the main runs: E. g., 
in simulation No. 7 about 20% of the vortices switched over to the other polarization, while in No. 6 no vortex did 
so (Fig. 3) . A closer inspection of the initial spin configurations shows that r depends both on the position of the 
vortex center within a lattice cell and on the dynamical state of the vortex. 

An additional test of the above conclusion was made by leaving out the first 500 time units of each main run 
(dash-dotted line in Fig. 3). Then we expect that the vortices have no memory of their initial condition and the 
resulting r should be the same as in the type-I simulations (within the statistical errors). In fact, this is confirmed 
by comparing No. 8 with No. 2 in Table I. 

V. CONCLUSION 

In this work we used a very simplified Hamiltonian for the cores of both planar and nonplanar vortices. Adding white 
noise to the local fields in which the classical spins precess we obtained a Landau-Lifshitz equation with multiplicative 
stochastic forces. The stationary solution of the corresponding Fokkcr-Planck equation exhibits two maxima for the 
two possible polarizations of the nonplanar vortex and a saddle point for the planar vortex, if the anisotropy parameter 
lies in a certain, temperature-dependent range. 

We calculated the rate n for the transition from one polarization to the opposite one. Our results were tested by 
long-time Langevin dynamics simulations of the full many-spin model at three temperatures well below the Kosterlitz- 
Thouless phase transition temperature. The agreement is rather good, considering that the vortex core was described 
only approximately by using only the four innermost spins. We did not make any tests for higher temperatures 
because the probability for the spontaneous appearance of a vortex-antivortex pair in the vicinity becomes too large; 
the interaction with this pair could then influence k. 

We emphasize that the above results were obtained by effectively averaging over many initial conditions. This is 
necessary because our simulations demonstrate that the transition rate depends very strongly on the initial condition, 
i. e. both on the position of the vortex center within a lattice cell and on the velocity of the vortex at this position. 
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Fig. 1. (a) Stationary probability distribution (14) (without normalization), using the simplified core Hamiltonian 
(26). TOi and ipi are the deviations (20) and (19) of the out-of-plane components and in-plane angles from their static 
values, resp.. 
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Fig. 2. A vs. /3 phase diagram. Above the bifurcation curve A(/3) there are two maxima in Vsti corresponding to 
the two polarizations of a nonplanar vortex, and one saddle point corresponding to a planar vortex structure. 
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Fig. 3. Percentage of vortices which have not yet made a transition to the opposite polarization up to the time t 
of a type-II simulation. The temperature is T = 0.15. The solid and dashed lines represent two sets of runs (No. 7 
and 6 in Tabic I) with different initial configurations (which arise from using two different lengths to for the pre- run) . 
The dash-dotted line results from sampling the transition times from the simulations No. 4 - 7 of Tabic 1, omitting 
the first 500 time units of each main run. 



TABLE I. Transition times r from simulations, with statistical errors Trms/T, compared to the theoretical estimates nh 



No. 


type 


T 


to 


t 


A^o 


N 


r 


Trms/r 


Tth 


1 


I 


0.1 


1200 


3800 


497 


477 


92516 


22% 


70334 


2 


I 


0.15 


1200 


3800 


407 


158 


4016 


6% 


3386 


3 


I 


0.2 


1200 


3800 


100 


1 


825 


10% 


743 


4 


II 


0.15 


1200 


4000 


254 


100 


4291 


8% 


- 


5 


II 


0.15 


2200 


4000 


264 


100 


4120 


8% 


- 


6 


II 


0.15 


3200 


4000 


181 


100 


6741 


11% 


- 


7 


II 


0.15 


4200 


4000 


405 


100 


2859 


6% 


- 


8 


II 


0.15 


- 


3500 


905 


400 


4286 


5% 


3386 
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